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TECHNICAL MEMORANDUM 1281 


UNSTABLE CAPILLARY WAVES ON SURFACE OF SEPARATION 
OF WO VISCOUS FLUIDS* 

By V. A. Borodin and Y. F. Dityakin 

The study of the breakup of a liquid jet moving in another medium, 
for example, a jet of fuel from a nozzle, shows that for sufficiently 
large outflow velocities the jet breaks up into a certain number of 
drops of different diameters. At still larger outflow velocities, the 
continuous part of the jet practically vanishes and the jet immediately 
breaks up at the nozzle into a large number of droplets of varying 
diameters (the case of "atomization") . The breakup mechanism in this 
case has a very complicated character and is quite irregular, with the 
droplets near the nozzle forming a divergent cone. 

Rayleigh (reference 1) was the first to make a theoretical study 
of the jet and to establish the possibility of droplet formation. The 
disturbance of a jet of an ideal fluid flowing into a vacuum and having 
a wave length 4.4 times as large as the diameter of the jet is shown 
to grow more rapidly than other disturbances; eventually, the jet 
breaks up into droplets of the same diameter. Rayleigh succeeded in 
determining theoretically the drop diameter, the value of which agrees 
well with tests on jets issuing with very small velocities. Later, 
the viscosity of the jet was also taken into consideration. The 
viscosity is found to decrease the rate of amplitude increase of the 
disturbances but the ratio of the optimal length of the wave to the 
diameter of the jet remains unchanged. 

Other authors that studied the conditions of the axial-symmetrical 
breakup of a jet of a viscous liquid found that the ratio of the optimal 
wave length to the jet diameter was somewhat greater than that computed 
by Rayleigh. 

In addition to the viscosity, Tomotika (reference 2) took into 
account the density and viscosity of the medium surrounding the jet 
and obtained good agreement with tests on jets issuing with very small 
velocities for which droplets of the same diameter are formed. 
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Neither of the aforementioned theories of the breakup of a liquid 
jet provided a basis for the phenomenon for the case of breakup into 
droplets of different diameter, a fact that is explained by the 
idealized conditions of the problem. This idealization consisted either 
in neglecting the viscosity of the jet, the density, and viscosity of 
the surrounding medium, or the inertial forces. Such simplifications 
were assumed in view of the complicated mathematical equation (generally 
transcendental) that determines the relation between the wavelength 
and the increment of the vibration amplitude. 

In the present paper, an attempt is made to provide a mathematical 
basis for the possibility of the appearance of droplets of different 
diameters as a result of the jet breakup on the basis of the considera- 
tion of unstable capillary waves on the surface of separation of two 
viscous liquids. 

For simplif ication of the solution of the problem, particularly 
for obtaining the algebraic characteristic of the equation, the 
lengths of the capillary waves on the surface of the liquid jet are 
assumed to be so small in comparison with the jet radius that the jet 
may be considered infinitely large; study of the stability of the 
plane surface of separation of two infinitely extending viscous fluids 
can thus be made. This assumption represents a considerable degree of 
idealization but nevertheless permits a qualitative explanation of 
not one but several unstable capillary waves that, in passing i through 
the jet, lead to the formation of droplets of differing diameters. 

The existence of several unstable capillary waves is demonstrated 
that can lead to the breakaway of several infinitely long strings of 
different dimensions from the partition surface. The problem investi- 
gated gives a rough approximation of the disintegration pattern of a 
liquid jet in another medium and- does not pretend to explain the com- 
plicated mechanism of the limiting form of the disintegration of a 
jet, namely, atomization. Nevertheless, one of the peculiarities of 
atomization, the appearance of a dimension spectrum of the droplets, 
begins to appear even for the given idealized consideration of the 
stability of the partition surface. 

1. Equations of small waves and their solution. - A plane surface 
of separation of two infinitely extending viscous fluids (fig. l) is 
considered. The viscosity and density of the lower fluid are denoted 
by p;l and P]_, respectively, and of the upper fluid by p 2 and P 2 . 

The lower fluid is assumed to move with the velocity and the upper 

fluid with the velocity V 2 , the direction of motion being the same and 
the velocities Independent of y. 
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A study of the character of the equilibrium of the surface of sep- 
aration under the action of the viscous forces and the forces of surface 
tension that impart to both liquids small disturbances parallel to the. . .. 
x-axis is presented. The fluids shall' be considered incompressible and 
weightless and shall cause certain disturbances to the components of the 
motion. 


V x - V + v x 


V = v 

y y 


p = p + p* 


It is further assumed that the velocities of the imposed disturb- 
ances and their derivatives up to the third inclusive are small and that 
the magnitudes of the second- order smallness may be neglected. 

From the Wavier-Stokes equations, the following equations of the 
imposed disturbances are obtained: 


( 1 . 1 ) 


Sv x ^ v 5v x 1 „ A 

p ax x 


c>v dv y 

__A. + y — -*■ 

at ax 


i a^ 

P Sy 


uA 


iV, 


where u = p/p is the kinematic viscosity. 
The equation of continuity is 

av Y a Vv 

5 + — L = o 

ax ay 


( 1 . 2 ) 


By introducing the stream function of the disturbance 




(1.3) 


and by eliminating the pressure p* from equations (l.l), the idealized 
equation is thus obtained in the Helmholtz form 


V AAty 


Y - 0 

ax at 


(1.4) 
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Let the stream function of the imposed disturbance be a periodic 
function of x and of the time t: 

* - (a = ^ (1.5) 

where a is the propagated circular frequency of the vibrations (the 
wave number)^ \ is the wavelength of the imposed disturbance, 

P = P r + iPj is the complex frequency of vibrations in time, P r 

is the real frequency of vibration in time, and p^ is the increment 
of the growth of vibration or the decrement of damping. 

The character of the wave motion on the surface of separation 
after the imparting of disturbances to both surfaces will thus depend 
on the sign of the imaginary part of the frequency p^. If p.j_ is 

positive, there will be an Increase in the wave amplitude with time; 
if Pj_ is negative, there will be a damping of the wave amplitude; 
finally, if P r = 0, there will be an aperiodic increase (p^ > 0) 

or a damping (P^ < 0) of the wave amplitude. By substituting expres- 
sion (1.5) in equation (1.4), the following equation is obtained: 

Uf IV - (2a 2 u - ip) f" - (iPa 2 - ua 4 ) f - iVa (f" - a 2 f) = 0 (1.6) 

The problem of the characteristic values of a homogeneous system 
of equations of the fourth order will be considered. 

p 

By setting f" - a^f = 'p , a system of equations of the second 
order is obtained. 

cp " + _ a 2 ^cp =0 f" - a 2 f = q> (1.7) 

Hereinafter, the following notations are introduced: 

/ P - V 2 a ~ 

= U 1 a -™ 2 ( 1 . 8 ) 


The solution of the first of equations (1.7) has the form 

cp = C 1 e 4m l^ r + C 2 e" in, iy (1.9) 

By substituting expression (1.9) in the second of equations (1.7), 
a non-homogeneous equation is obtained for which the solution is 
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e irQ iy 

f p p Cl - 

ran * + op 


C 2 + e a y C 3 + e" a y C 4 

m-, + ar 


(1.10) 


The stream function for the lower and upper liquids according to 
equation (1.5) will be 


V i = e 


i(cuc-Pt.) 


e lm iy 


2 2 1 

m-, ^ + a' 1 - 1 - 


C, - 


.- lm iy 


m^ + a‘ 


2 C 2 +e<W C 3 + e “ ayC 4 


(1.11) 


y = e i (ax-Pt) 


>im2y 


mo + ar 


c 5 - a 6 ♦ .<* c 7 + e-“ y c 8 ) (1.12) 

m 2 c + a 


The arbitrary constants C 4 must be determined from the conditions 
on the surface of separation and at infinity. 


2. Boundary conditions. - The boundary conditions of the problem 
will be as follows: 


1. At infinity (y = +“0, finite solutions must be maintained for 
and ^ 2 . Hence, the arbitrary constants of the terms with positive 

exponents for V 1 and with negative exponents for V 2 must be equated 
to zero: C 4 = C 3 = Cg = Cg = 0. Thus, equations (1.10) and (l.ll) will 
have the form 


y = e i(ax-pt) 


y 2 = (a-x-Pt ) 


e ~ tm iy 

m 4 2 + a 2 



e im 2y 

2 , _2 
m 2 + a 




( 2 . 1 ) 


2. On the surface of separation, there must be no slip, that is, 


M y.O ‘ My,0 

( v ,yi) y .o = ('Vs) 7= o 


or 
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f*i\ 

V oy / y=:0 V Sy/ y= 

f fjM = f ^ 

V Sx/ n V ox / V=| 


y=o J 


( 2 . 2 ) 


3. The tangential stresses on the surface of separation are 
continuous 


^l< A Vl)y=0 = p 2< A V 2 )y=0 (2.3) 

4. The difference between the normal stresses p y ]_ and p y2 

on the surface of separation is equal to the pressure brought about by 
the surface tension; that is. 




p y2 




Pi + 2p 1 




p 2 + 2p 2 


Sv 


y2 


ay 


\ S 2 h 

)= " T ~ 

/ ox 


(2.4) 


where T is the capillary constant of one liquid relative to the other 
and h is the rise in the surface of separation at the point x. 

By using equations (2.1), the boundary conditions (2.2) are obtained 
in the form 


im]_ 


+ a‘ 


- C 2 - aC 4 + 


'2 "5 

2 „2 4 + ™ 2 , „2 


im 2 

— g 2 C 5 - aC 7 = 0 

m 2 + a 

C 


- C 7 = 0 


(2.5) 


m]_ + a + a 

Similarly, the boundary condition (2.3) is obtained in the form 


Pi 


/ c N 

l 2 . “l 2 „ _ 2 „ 

p 2 


( _ 2 , _2 
V + cc j 

1 Q/ * p p ^ '-M 

' m n c + ar 


+ C 7 ]a 2 + 


rtin ^ + a ? 


Dip p 

Z C q + arC- 

2 2 5 1 
mo + or 


( 2 . 6 ) 
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The pressures p-j_ 
(2.1). Thus 

and p 2 

are computed 

from equations (l.l) 

Pl - P 1 e 1(t “' Pt) 

"oe-imiy i 

r v ! a 2 

iV x a 

. jjP -f TTli Ui I Co + 

m]_ 2 + a 2 ' 

^ m l 

mi 

mi x x ) 


(2.7) 


p 2 = pgeK^-Pt)! 


+ il^a 2 - V 1 a - C^J 

^ + IE ♦ As ♦ (» - V 2 a) ^ C 7 


[_m 2 2 + a 2 \ m 2 


mg mg 
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The rise of a point on the surface is a periodic function of x 
and t . 


h , He 11031 - 1311 ) (2.8) 

where H is the maximal rise of a point on the surface of separation. 

The velocity of the raised point on the surface of separation is 


(V) y , 0 'fe) 0 -| 

y=0 


Sh + y Sh 
ot dx 


(2.9) 


After differentiating expressions (2.1) and (2.8) and by substi- 
tuting in expression (2.9), the following equation is obtained: 


H = — - [ C 4 - 2 


aY 1 - (3 


m-i 2 + a 2 


( 2 . 10 ) 


By substituting equation (2.10) in (2.9) and by differentiating 
equation (2.9), 


o^h 


5x 2 = ^1 " P U 2 + a 2 4 


- C,)e i ( ax -P t ) 


(2.11) 
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By computing the derivatives dv^/Sy and Sv^/Sy and substi- 
tuting them simultaneously with expressions (2.7) and (2.11) in (2.4), 
the following boundary condition is obtained: 


a 


m-L + or 


-hl a2 + iP]V]_a - iPp-^ 


m n 


mini 


Tar 


aV 1 - p 


p-j_3 - Pi^a, - i2[ijar + 


aV n - 3, 


C/i + 


a 


mo + ar 


P2°- 2 - iPgVga. + iPPg 


3m 2P2 


C^ + (pg3 “ P2^2 a — 0 


The following nondimens ional parameters are then introduced : 


U 1 a 




E 2 = TTE 
2 


A = 


V. 


K = 


2 

tp x 

2~ 

Ul a 



( 2 . 12 ) 


(2.13) 



NACA TM 1281 


9 


where c = p/a is the complex wave velocity. 

Equations (1.8) can then he represented in the forms 


m-j_ = a i(Z - E^) - 1 


mg = a t\j i(ZA - Eg) ~ 1 


Equations (2.5), (2.6), and (2.12) are represented in nondimen- 
sional parameters. The following notations are first introduced: 


a l* " ^1 Z-B n 


2 N 

/\/i (Z-R-l) - 1 + Z " R 1 


= PtS 


1“1 


V - ( Z - E r 21 + s^z) * 


Vi 


2(1 - 2Rg + 2AZ) 2 

Cl ~ (ZA-Rg)^/l(ZA-Rg) - 1 ^ 2a Cl 


a 3* = " — 1 


a 2 Z " R 1 a 2 


d x * = p 2 a 2 (ZA-R 2 + 2i) = P 2 a- 2 d 1 


1 i 

° 3 * * ' ~3- ^ 


>* " I VUZ-HP - 1 = ^ 

(ZA~Rg) - 1 Cg 


Co* = 

2 a ZA-Bg 


a 


(2.14) 


where a-j_, ag, a 3 , h^, C]_, Cg, c 3 , and d^ are likewise nondimen- 
sional magnitudes. 
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The following system of equations is then obtained for the con- 
stants C 2 , C 4 , C5, and C 7 : 


a l*C 2 


b 1 *C 4 + c 1 *C 5 + d 1 *C 7 = 0 


a 2*^2 " “£4 + c 2*^5 ~ a ^ 7=0 

-a 3* C 2 + c 4 + c 3* C 5 ~ C 7 = 0 


> 


(2.15) 


K C 2 - C 5 = 0 J 

This system of homogeneous equations has solutions different 
from zero if its determinant is equal to zero. By setting up the 
determinant and expanding 


2K(a-^ + c-^) + (d^ — Kb-^) (a 2 + iCc 2 ) + (Kb^ + d-^) ( Kc ^ — a-^) = ^ 


By solving this equation for Z, the following wave equation of 
the 18th degree with complex coefficients is obtained : 

r 18 z18 + ( r i7 + is i7) z17 + ••• + ( r 1 + is x ) Z + (r Q + is 0 ) = 0 

( 2 . 16 ) 

The real and imaginary parts of the coefficients depend on the 
five nondimens ional parameters: R-j_, R 2 , A , R, and K. 


5. Investigation of roots of characteristic equation, - The 
increase in oscillation, that is, the loss of stability of the sur- 
face of separation, arises from those waves for which the imaginary 
part of the frequency is positive (P 4 > 0). Hence, the investiga- 
tion of the roots of equation (2.16) should determine those ranges 
of the parameter N or the wave number a in which the complex 
roots of the equation lie in the upper half -plane. 

By the Rayleigh hypothesis, the further development of an 
unstable deformation, that is, the form and dimensions of the parts 
breaking away, is determined by the critical (or optimal) disturb- 
ances. The critical disturbances may be defined as those that 
develop more rapidly than the others or that correspond to the 
maximum, increment of the growth Pi. This principle of deter- 
mining the character of the unstable deformations by the character 
of the maximum unstable disturbance has been experimentally 
confirmed by a number of investigators (reference 3). 
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In the case considered, the growth in the amplitudes of the- 
oscillations will lead to "breakaway of infinitely long strings 
from the surface of separation, similar to the formation and • 
breakaway of wave crests. The separation will take place for such 
values of a or wavelengths \ for which 8^ has the maximal 
value . 

If a spectrum of small-period disturbances that can be developed 
into a Fourier series can be assumed to be imposed on both liquids, 
the harmonics with the wavelengths equal to the wavelengths of 
the maximal unstable disturbances bring about a separation of 
infinitely long strings from the partition surface. Because the 
characteristic dimension (for example, the diameter of the transverse 
string) is connected with the length of maximal unstable disturbance, 
strings of different dimensions will break away from the surface of 
separation. In figure 2, the scheme of formation of such strings 
for three successive instants of time is shown. 

Investigation of the roots of the simplest particular case of 
equation (2.16) is presented. 

Let both fluids be stationary and their kinetic viscosities 
the same. In this case, V-j_ = V 2 = 0, b-j. = u 2> m l = m 2> A = 1, 

% ~ r 2 ~ an< ^ equation (2.16) goes over into an equation of the 

8th degree whose coefficients depend only on the two parameters 
K and W : 


A q Z 8 + (A ± + iB 1 )Z 7 + (Ag + iB 2 )Z 6 + (A^ + iB 3 )Z 5 + (A^ + LB 4 )Z 4 + 

(A 5 + iB 5 )Z 3 + (Ag + iBg)Z 2 + (A ? + 1B ? )Z + Ag = 0 (3.1) 
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where 

A 0 = - (1 - K) 2 

A 2 = 2K (K - 1) N - K 4 + 2K 3 - 4K 2 + 6K + 13 
A x = - 2K (K - l) 2 
A 3 = 4K 2 (K - 1) N - 2K (3E 2 + 13) 

A ? = 2K 3 N 2 

A 4 = - K 2 N 2 + 2 (E 4 - K 3 + 3K 2 + 5K) If - 12K 3 + 26K 2 - 10K - S 

A 5 = - 2K% 2 + 12K 3 If - 8E 3 - 8K 2 

A e = - K 2 (1 + 2K) If 2 

Ag = (1 - K 2 ) K 2 N 2 + (4K 4 - 10K 3 - 4K 2 - 6K) N - 8K 3 + 12K 2 - 8K + 4 

B x = 2 (K 2 + 2K - 3) 

B 3 = 3E 4 - 14K 3 + 13K 2 + 18K + 13 - 81® 

B 2 = 4K (K - l) 2 
B 4 = QK 3 - 4IC 2 - 2 OK - 2K 3 N 
Bg = 4K 2 (1 - K) If 

B 5 = - 2K% 2 + [2K (1 + K) (1 + K - K 2 ) + 8K 3 + 4K 2 + 4K]lf + 

4(K - 1 - K 2 ) (1 + If - K 2 ) + 4K 4 - 2 OK 3 - 4(K - 1 - K 2 ) 2 

B ? = [K 2 (1 + K) 2 - 2K 4 ] If 2 + [4K 4 + 4K (1 + K) (K - 1 - K 2 )] If 

(3.2) 
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The characteristic equation (3.1) is a polynomial whose 
coefficients depend nonlinearly on the two parameters K and N. 

Each pair of values of the. parameters K and N or each point 
of the. plane KN correspond to the completely defined polynomial 
(3.1) , that is, completely determined values of the eight roots 
of the polynomial. In the plane KN, it is evidently possible to 
find a curve, each point of which corresponds to the polynomial 
(3.1), that has at least one root located on the real axis so that 
only in crossing this curve is a crossing of the roots through the 
real axis possible. This curve breaks up the plane KN into 
regions, the points of which each correspond to polynomials (3.1), 
that have the same number of roots with positive imaginary part. 

These curves are constructed by making use of the method of 
Y. I. Neimark (reference 4) that permits a breakup of the plane 
of the parameters for the roots of the polynomial lying in the left 
or right ha If -plane. 

The substitution Z = -it, is made. The upper half-plane of 
the roots of equation (3.1) is transformed into the left half -plane 
of the roots of the equation 

- A 0 ^ + i (A-j_ + iB^) t? + (Ag + iBg)£ 8 - i(A 3 + iB 3 )£ 3 ~ (A 4 + 1 ^ 4 )^ + 
i(A 5 + iB 5 R 3 + (Ag + IB 6 )C 2 - (A ? + iB ? H + A 8 = 0 

(3.3) 

By substituting £ = i£/t| in the preceding equation and multi- 
plying the result by tj 8, equation (3.3) is reduced to the form. 

F(£,ti) + iG(S,n) = 0 (3.4) 


where 


F(£,ti) = A 0 £ e + A l £ 7 t| + AgfcV + A 3 £ 5 T ] 3 + A 4 £ 4 t ) 4 + A 5 £ 3 t ) 5 + 

A 6 £ 2 q 6 + A 7 ir \ 7 + A 8 q 8 (3.5) 

G(£,t,) = B^q + BgfcV + B 3 ^T1 3 + B 4 S 4 B 4 + B 5 I 3 ti 5 + B 6 e 2 B 6 + B 7 £q 7 



miiiinii i iiiiiii 


ii ii mil mi i inn n i 


I II INI 


mu 1 1 i in i mi i 


1 1 ii i 
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If R 2n ia s P ace °f complex polynomials of degree n and 
D(k,n - k) is the manifold of polynomials R 2n having k roots to 

the left and n - k roots to the right of the imaginary axis of the 
complex sphere, then by setting up the following table: 



(3.6) 


and by making the transformation 

Aq + ^1 B 1 A 1 + A 2 + A 3 + X 1 B 4 ...A f 




B 2 


. . .0 


(3.7) 


table (3.7) is found to correspond to a polynomial of the same type 
with respect to the distribution of the roots relative to the imaginary 
axis, as in equation (3.4). 

From table (3.6), an inequality . is obtained that defines the 
region in the plane KN corresponding to the presence of the first 
root of equation (3.1) in the upper half -plane: 


A 0 B 1 < 0 (3.8) 

By setting = - Aq/b 4 in table (3.7) 

( A i B i - a 0 b 2 )/b 1 ( a 2 b 1 - a 0 b 3 )/ b 1 ( a 3 b 1 - A q b 4 ^ / B i ••• A 7 A el 

B 1 B 2 B 3 • • - b 7 0 

(3.9) 


Because A-jB^ - AqB 2 = - 16K(K - l) 3 < 0 for K>1, by multi- 
plying the elements of the first rows of (3.9) by B- l 2 /( a 1 B- l - AqB 2 ) 
and changing signs in the second row 


B 1 B i B 2 D 3 
- B 1 - B 2 - B 3 - B 4 





> 


(3.10) 
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where 


Wi 

•Vs) 

A 1 B 1 - 

a 0 B 2 

B 1 ( A 3 B 1 

" AqBa) 

A 1 B 1 " 

A 0 B 2 

( a 4-Bi 

- AqP 5 ) 

A 1 B 1 - 

A 0 B 2 

B l( A 5 B l 

“ A 0 B 6^ 

AiBi - 

a 0 B 2 

B l'Vl 

- Vt> 

A 1 B 1 - 

a 0 B 2 


Dl 


= D 2 


= L, 


= D/i 


= D c 


B 1 A 7 


A 1 B 1 “ A 0 B 2 
B 1 A 8 

A 1 B 1 “ A 0 B 2 


= D c 


= D- 


(3.11) 


The first row of table (3.10) is left unchanged hut to the 
second row is added the first row. Thus 

B i D 1 d 2 D 3 D 4 D 5 B 6 B 7 

0 D i _ b 2 D 2 “ b 3 D 3 “ B 4 D 4 “ B 5 D 5 ” B 6 B 6 “ B 7 C 7 

(3.12) 

From the preceding calculations, an inequality is obtained that 
defines the region in the plane KN that corresponds to the presence 
of the second root of equation (3.1) in the upper half -plane 


B 1 (D 1 - B 2 )< 0 


(3.13) 
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By carrying out a ■ transformation, similar to (3.7) of table (3.12), 
B-^ + X 2 (®l “ Bg-) ^ Xg (Dg ” B 3 ) Bg ^2 ( B 3 "" B 4)* ,,B 7l 


P 1 " B 2 


P 2 " B 3 


..C- 


► (3.14) 


By setting Xg = - B l/(B-] ~ Bg) and substituting in (3.14) 
D 1 (D 1 - Bg) “ B q (Bg - Bg) Bg ( P l ” ®2 ) “ B q( P 3 “ B 4 ) 


D 1 “ B 2 
P 1 " B 2 


P 1 “ b 2 
d 2 - b 3 


D 3 (Di - Bg) - B 1 (D 4 - B 5 ) 


K (3.15) 


P 1 “ b 2 
d 3 - b 4 


Because D 1 (D 1 - B 2 ) - B 4 (D 2 - B 3 ) > 0 for K>1, by multi- 
plying the elements of the first row of (3.15) by (B-^ - Bg) 2 / 

[%(% - B 2 ) - B l( D 2 “ B 3)] 


Di - B 2 


P 1 “ B 2 


[D 2 (D X - Bg) - B^Ds - B 4 )](D x - Bg) 
Dp (D x - Bg) - B 1 (D 2 - B 3 ) 


(3.16) 


b 2 - b 3 


The elements of the first row are subtracted from the elements 
of the second row of table (3.16). 


^j~D 4 - Bg £(Bg < P 1 ” B 2^ “ B ].( P 3 ” b 4)J ( P 1 “ b 2^ • • •""! 

S 0 (D - B ) - [ p 2( p l - B 2) - B l( p 5 - B 43 ( D 1 - B 2 ) | 

2 3 D 1 (D 1 - Bg) - B 1 (Dg - B 3 ) ***- 

(3.17) 

From the preceding table, an inequality is obtained that defines 
the region in the plane of the parameters KN that corresponds to the 
presence of the third root of equation (3.1) in the upper half -plane. 


(Dp - Bg) 


* > [P 2 ( p i - Bg) - B^Dg - B 4 )] (£1 - Bg) 
( 2 “ 3 D 1 (D 1 - B 2 ) - B 1 (D 2 - B 3 ) 


< 0 

(3.18) 
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Similar conditions can be obtained for all the remaining roots 
of equation (3.1). This investigation has been limited to the three 
conditions that are sufficient for proving the existence of several 
unstable waves. 


By replacing inequalities (3.8) , (3.13) , and (3.18) by equations, 
the equations of the curves determining the breakup of the KN plane 
into regions are obtained. The most interesting case of large 
K = i^/Vg^l is considered. From inequalities (3.8), (3.13), 


and (3.18) and by considering equations (3.2) and (3.11) and neg- 
lecting small powers of K, the following equations are obtained : 


2 (K - ir(K + 3) = 0 
CqN 3 + e^N 2 + e^N + 63 — 0 


4(K 2 - 1)(K + 3)N + K(K - l) (K 3 + 17K 2 - 96K + 99) = 0 


where 


(3.19) 


e Q = 128 (K 4 - K 3 - 23K 2 - 39K - 18) 


= 592K(K 5 + 8.4K 4 + 3.18K 3 - 96K 2 - 20. 3K + 0.98) 


e 2 = 9K 2 (K 6 + 8.4K 5 - 97. 3K 4 - 2045K 3 + 1700K 2 + 390K + 363) 


e 3 = 24K 5 (K 5 + 12. 3K 4 + 306K 3 - 4100K 2 + 12,300K - 7000) 


By plotting the curves (3.19) in the KN plane and separating 
by hatched lines the regions corresponding to the signs of the 
inequalities (3.8), (3.13), and (3.18), the diagram shown in 
figure 3 is obtained. This diagram shows that for K> 0 and N>0 
a region of values of K and N exists that corresponds to the 
presence of three roots with positive imaginary part, that is, 
of three unstable waves on the surface of separation. 


The division of the KN plane for the remaining roots could 
establish regions with a still greater number of roots with posi- 
tive imaginary part. The given Incomplete diagram already shows, 
however, the existence of several unstable waves. In the presence 
of a maximum (3^ or c^, several infinitely long strings will 
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■break away from the surface of separation, the cross-sectional 
dimensions of which will depend on the wavelength of the critical 
disturbance. 


Translated by S. Reiss, 
National Advisory Committee 
for Aeronautics. 


REFERENCES 

1. Rayleigh: The Theory of Sound. Dover Pub., 2d ed., 1S45. 

2. Tomotika, S.: On the Instability of a Cylindrical Thread of a 

Viscous Liquid Surrounded by Another Viscous Fluid. Proc. 

Roy. Soc. London, vol. CL, no. A870, ser. A, June 1935, 
pp. 322-337. 

3. Petrov, G. I.: On the Stability of Turbulent Layers. Rep. 

No. 304, CAHI, 1937. 

4. Neimark, Y. I.: On the Problem of the Distribution of the Roots 

of Polynomials. DAN, T. 58, No. 3, 1947. 



NACA TM 1281 




3 


Figure 2. 


K 


K 


I 


Figure 3 


NACA-Langley - 4-18-51 - 075 


NASA Technical Ut 


3 1176 01441 1822 






